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ABSTRACT 

It is noted that a surface layer of matter in contact with a sufficiently super-Eddington, radially combed 
photon flux typically attains a relativistic coasting state whereby the radiation does not accelerate the 
matter. The final bulk Lorentz factor of this matter is therefore determined by geometry. Radiation 
that scatters off this layer is most likely to be observed along the velocity vector of the matter, where it 
' would be most strongly polarized. 

■ Subject headings: black hole physics — gamma-rays: bursts and theory 

(N ; 

A remarkable feature of cosmologically distant gamma ray bursts is their high isotropic equivalent luminosities, which 
can exceed the Eddington limit by up to 14 orders of magnitude. Such a flux would sharply accelerate baryonic matter 
, that stood in its way and, even for a small quantity of such matter, most of the GRB energy would be expended on its 
' acceleration. This would require salvage mechanisms such as internal shocks, which, in order to be effective in converting 
\l ' efficiently to gamma radiation, must occur near or beyond the photosphere, i.e. at large radii. 

Is it necessary to have any baryonic matter at all in GRB? If GRB emerge along field lines that thread an event horizon, 
they need not receive any baryons, and most of the pairs could annihilate at relatively small radii when the temperature 
falls below 2 x lO^K. The indirect evidence for matter in the fireballs is a) non-thermal spectra, which suggest that much 
of the energy passes through charged particles at low optical depth, and b) the existence of afterglows, which would not 
be generated purely by gamma rays. There are, however, models that could yield non-thermal emission and afterglow 
, without a significant amount of baryons. For example, fireball energy is easily converted to ultrahigh-energy particles 
' if neutrons decay within the baryon-poor jet after creeping in across field lines (Eichler & Levinson 1999, Levinson & 
' Eichler 2003). The ultra- high energy particles could carry enough energy per electron that they could deliver the required 
non-thermal gamma ray energy and afterglow energy while remaining optically thin to the non-thermal gamma rays. 
' Other models posit pair creation purely from gamma rays (Thompson & Madau 2000, Beloborodov 2001) or Poynting 
(— I . flux (Lyutikov, Parlev, & Blandford 2003) at large radii. The photosphere, which would be at large radius for a baryon 
^ i' dominated plasma, can in view of these non-baryonic alternatives be at much smaller radius. To quantify this, note that 
' I if the energy that remains kinetic - ultimately destined to power the afterglow - is carried as normal electron-ion plasma, 
^ . then the photosphere occurs at r = 1/(1 — /3)crT'^e 10^^-^b5o( i^)~'^(o^)^^ cm. Here L^so is the kinetic luminosity of the 
' baryons, ATrT^pmprieC^d^r^ , in units of 10^'^crg/s, the neutron to proton ratio in the ions is assumed to be unity, and we 
^ , have tentatively included the possibility that the internal Lorentz factor for the ions may be larger than unity. These 
' considerations would together suggest that ri2 > 1 and that F > 10^, which is widely though not universally believed for 
• I GRB fireballs. If, on the other hand, the fireball plasma is mostly pairs then the photosphere can be dictated by pair 

' recombination and occur at considerably smaller radius, {vph < lO-'^^cm), e.g. Eichler (1994), Eichler & Levinson (2000). 
'H ' A now familiar suggestion for getting kinetic fireball energy into gamma radiation, is via internal shocks near (Eichler 
■ - - ' 1994) or downstream of (Rees & Mezsaros 1994) the gamma-ray photosphere. The question still remains whether the soft 
gamma ray emission, is emitted from optically thin regions or from a photosphere, or both. The low frequency parts of 
GRB spectra are frequently too steeply rising, it is reported (e.g. Preece et al. 2002), to be consistent with optically thin, 
efficient synchrotron radiation. Alternative pictures, which might allow a steeper rise in the soft part of the spectrum, 
involve Compton upscattering (Shaviv & Dar 1995, Lazzati et al. 2000) of softer photons by the energy-bearing, relativistic 
outflow. Here too, however, the bulk and internal Lorentz factors must be fine tuned to get the right peak energy. 

However, there is another alternative: that radiation from or near a photosphere is the primary energy source, this 
photosphere is largely determined by net annihilation of positrons - which which loosely associates the peak energy with 
nieC^ - and baryon involvement, while generating additional effects, occurs too far downstream to spoil this association. 
Recent evidence of GRB collimation suggests that the fireball, though sufficiently baryon-pure to allow high bulk Lorentz 
factors, is surrounded by baryon-rich material either the envelope of a host star, or a wind from the accretion disk 
(Levinson & Eichler 2000). This allows an entirely different set of models for GRB's, one in which the GRB fireball 
would be pure photons (or Poynting flux) except for matter that is "shaved off the inner wall of the baryonic sheath 
that coUimates the fireball. The shavings could be accelerated to relativistic Lorentz factors, and could play whatever 
role is ascribed to matter in GRB fireballs, such as eventually supplying non-thermal particles and generating afterglow. 
However, they need not be the primary source of soft gamma rays and they need not be subject to constraints on F 
derived from time variability. 

The purpose of this Letter is to consider the interaction between a radiative fireball and a sheath of baryonic matter 
when the latter coUimates the former. It is argued that the moving wall of matter that is in contact with the fireball 
maintains a Lorentz factor F that is the inverse of the sine of the impact angle between the photons and the outflow. This 
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is due to the fact that this matter "sails" on the radiation pressure. This picture could yield an estimate of the Lorentz 
factor of the matter that is associated with GRB fireballs. It also predicts that the polarization is typically high; it can 
be close to 100 percent, and can exceed 60 percent with significant probability. This is considerably higher on average 
than in previous suggestions (e.g. Shaviv & Dar 1995, Eichler & Levinson 2003). 

Compton equilibrium between matter and radiation pressure has previously been considered in the context of an 
accretion disk by Sikora and Wilson (1981) Phinney (1982), Sikora et al (1996). Begelman and Sikora (1987) considered 
polarization from scattering off outflow, but did not specifically consider the correlation of such polarization with Compton 
equilibrium. Beloborodov (1998) considered the polarization vector in the case of a 1-D outflow from a disk with Compton 
equilibrium at modest /3, /3 ~ 0.5. The polarization he found is along the outflow axis, orthogonal to what is expected in 
the context of the present paper. The polarization effect proposed here is a consequence of high Lorentz factor combined 
with Compton equilibrium. 

Consider a flux of where I ~ 10^*. Suppose an electron ion pair is injected at radius Rq = XoRsch, where Rsch is 

the Schwarzchild radius. As shown below, the electron proton pair is accelerated radially to an asymptotic T of (Z/.Xo)^/^ 
when r 1. The fractional power 1/3 is due to the redshift suffered by the radiation in the frame of the accelerated 
matter. Note that for I ~ 10^^ and x ~ 10^, the asymptotic F is of order 200, quite consistent with popular estimates. 

At a collimating wall, however, the baryonic outflow can be constrained to be non-radial, and the radiation pressure 
can be more effective because there is less redshift. Though somewhat counterintuitive, this principle is similar to the 
fact that a ship sails faster when kept oblique to the wind by its rudder than when traveling directly downwind. If a 
collimating wall keeps the flow moving oblique to the radial direction (which we are assuming here is about the same as 
the direction of the photon wind) then the acceleration or deceleration can be far more effective. The extremely large 
super-Eddington factor 1 implies that the final Lorentz factor is determined by the coUimation geometry. 

The radiation force on an electron proton pair along its direction of motion /?, which moves at angle x relative to the 
radial direction f, is most easily calculated in the frame of the pair, where the cross section is the Thompson cross-section. 
In this frame the component of force F' along the direction of flow is given by 

F'-p = a-rT' ■ ^/c = (^T j d^'k' I'(z^', k')(k' • /c (1) 

where J^' is the comoving radiative flux, ut is the photon scattering cross section which throughout is assumed to be the 
Thomson scattering cross section. The acceleration along the wall in this frame, a', is a' = F' • /3/m. In the limit of strong 
radial combing of the radiation, /'(i^', k') can be factored into Tl{i'')S'^{k' — f'). Using the invariance of I' {u')dv' / v'^ , 
we can write 

-r-, ^ 5^{k-r) v'^ 

where f' is the k vector of the radially combed photons in the primed frame and Tvdv = ILedd/ ^T^r^ ■ In the lab 
frame, d(3T/dt = V^a = a'. The primed variables can be expressed in terms of the unprimed variables as follows: Using 
the aberration formula, 

cosx' = [cosx -/?]/[!- Pcosx] (3) 

and 

u' = uT{l-pcosx), (4) 

we can write 

drp/dt = ^Wt SHk-f) _ ^ )4 (ca.X -/?) 

'^^ Anr^mp cS^k' - f') ^ {1 - pcosx) 

When r 1 and 1 — cosx ^1-/3, the quantity cosx' = (1 ~ Pcosx)~^{cosx — /3) ~ 1, and S^{k — f)/S^{k' — r') oc 
[d{l - y ■ r']/[d{l - fc • r] = (1 + - /?), so 

drf3/dtr~.dT/dt = -{4T)-'^ld{l/x) (6) 
and the solution as a function of radius xRsch is 

r{x) = [-{l/xo-l/x)+Tl]'/'. (7) 

This solution could be applicable to loose material that suddenly finds itself in the path of the fireball, such as a decaying 
pickup neutron that has drifted in from the walls. For Xg ~ 10®, corresponding to a radius of ~ lO^^cm, and I ~ 10^^ 

to 10^**, this yields an asymptotic T of several hundred, which is an acceptable value, though Poynting pressure could 
drive it to an even higher value. Also, the Eddington luminosity increases with the number of pairs per baryon, so 1 could 
effectively be higher than 10^'^. 

However, for l/x ^ T'^, F grows very rapidly, possibly until the condition 1 — cosx ^ 1 — /? is no longer satisfied. If 
the acceleration were to lead to 1 — cosx 3> 1 — /3, it would follow that 1 — Pcosx 1 — cosx, and the deceleration would 
then be extremely powerful. For l/x 3> 1, the value of F is rapidly reduced to the value where cosx = /5j i-e. F = 1/sinx- 

It is of course also possible that the condition 1 — cosx "C 1 — /3 is never satisfied, e.g. (3 could be small and 1 — cosx 
could be of order unity. In this case the acceleration along the wall is even larger than given by equation (6) because 



3 



the radiation is not rcdshiftcd nearly as much; the Doppler factor r(l — f3cosx) that would have lowered the flux in the 
fluid frame when 1 — cosx is negligible is now less devastating. This only strengthens the conclusion that the material is 
accelerated until (and only until) (3 ~ cosx- 

The final value for V is thus 1/sinx- By the same argument, thermal motion is strongly Compton cooled. If the 
source has a finite size, then the incident photons have some finite angular spread in the frame of the fluid. However, any 
Compton heating this may cause is modest, because the frequency in the fluid frame v' is only i//r <C mgC^. Moreover, 
if the cylindrical radius of the fluid r is much larger than the source size Vg, as we have assumed here, then the angular 
spread of the photons in the fluid frame can be shown to be small, of order rg/r, (though larger than the angular extent 
of the source as seen from the point of photon-fluid impact). A more general treatment of flnite source size will appear in 
a separate paper (Levinson and Eichlcr 2004). 

The incident radiation thus makes an angle of 7r/2 relative to the direction of motion, as would wind relative to the 
motion of a frictionless sailboat whose motion is kept by its rudder oblique to the wind direction. For the extremely high 
super-Eddington factor 1 that is typical of GRB's, we believe that viscous drag by the inner layers does not change this 
conclusion significantly when the mean free scattering time of an ion is longer than the acceleration time to P — cosx- 

The radiation, when scattered off the matter, is beamed most strongly in the direction of motion, which corresponds 
to 90 degree scattering, and thus full polarization along the direction n x /3. This should be contrasted with the situation 
corresponding to the "head-on" approximation as termed by Begelman and Sikora (1987) in which the radiation in 
the frame of the scattcrer is moving almost parallel to the lab frame velocity as seen by the scattcrcr. The head-on 
approximation would be natural when the flow is ultra- relativistic and the angle of incidence is physically unrelated to T. 
In this circumstance, beaming along the direction of motion gives no polarization, and it is only the flnite spread of the 
beamed radiation that allows polarized contributions from scatterers moving slightly off the line of sight. For a hollow 
cone geometry, the maximum intensity is seen by observers with lines of sight on the cone. The maximum polarization 
they may observe is 25 percent. The polarization vector is perpendicular to the cone, i.e. each fluid element with direction 
f3 contributes polarization in the direction {h — (3) x n. Sufficiently far away from the cone, i.e. sufficiently far into its 
interior or exterior, the polarization can be larger than 25 percent, but only at the expense of lowered intensity. Thus 
the Vmax associated with large polarizations is relatively small and the probability of viewing polarization more than 75 
percent, say, is only 0.2 (Eichler and Levinson 2003). 

By contrast, scattering from a Compton sail reverses the roles of high and low polarization, so the largest Vmax is 
associated with high polarization. In the limit of a line of sight of 59 from the edge of an infinitely thin cone, where 
59 ^ and l/F 9o, most of the contribution to the line of sight comes from fluid elements with (3 near n, i.e. 
l3-nr^59. 

To set up a formal calculation, consider scattering off a thin cone of opening angle 9o defined by the azimuthal angle 
(p running from to 27r, and the line of sight is offset from the cone by angle A. Let us use coordinates such that 

the unit line of sight vector n = [sin{9Q + A),0, cos{9q + A)] is in x,z plane. The ingoing photon vector is ki = 
[sin{9o + Y)cos(j), sin{0o + Y)sin(l),cos{0o + ^)], the outgoing photon vector, k2 is parallel to n. The velocity vectors of 

the outflow, which form a cone around the z axis, are defined by /3 = [sin9QCos(l), sin9Qsin(f>, 0056*0]. 

In the fidly general 3-D scattering problem, the degree of polarization depends on the scattering angle w', where 
cosw' = k[ ■ whereas the Doppler boosting of the intensity of the scattered radiation depends on 9 = arcos{(3 ■ n) 
as I{u) = I'{h'')[T{l — l3cos9)]~''. Here k = 3 — a where a is the spectral index. The difficulty is that w' and 9' are 
defined relative to different axes. However, in the case that the incoming photon vector makes an angle of cos~^j3 with 
the velocity vector, the frequency in the frame of the scatterer is given by 

(8) 

and we assume that v' is conserved in the scattering. We can use the invariance of fci^fcg to write 

1 - cosw' = r^ — (1 - cosw;) = r2(l + (3cos9'){l - cosw). (9) 

1 + Pcos9' = 1 + P^^^^ = \ (10) 

1 — (jcosf) F^(l — PcosO) 

we can write cosw' in terms of laboratory coordinates as 

, , (1 — cosw) (cosw — /3cos9) 
(l-Pcos9) = (1-PCOS9) - 

The Doppler boost and polarization expressions from emission from a thin cone of opening angle 60 can now be expressed 
as functions of ^, F and 9o using the basic expressions 

cosw = sin{9o + A)sin{9o + ^)cos(/> -|- cos{9q + A)cos{9o + ^) (12) 

and 

Pcos9 = 0-h = l3[sin9osin{6o + A)cos(j) + cos6'ocos(6'o + A)] (13) 



Using the aberration formula, 
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In the limit of l/F <C and A <^ Og, we can expand the trigonometric formulae and obtain 

1 - cosw^ [{1/r - +4)'^]/2, (14) 

where cp = sinOocf), and 

l-/3cos6''-(l/r2 + A2 + <^2^/2 (15) 

whence 

cosw - pcosO ~ A/r (16) 

Thus, in the viewing direction A = 0, cosw' = and the radiation scattered into this direction after a single scattering 
is nearly fully polarized. It can be further shown in the single scattering limit that when F ^ 1, the polarization vectors 
of all the rays scattered into a line of sight with A = have polarization nearly perpendicular to the cone axis as seen 
projected on the sky, so that their sum is nearly fully polarized as well. 

In the limit l/T 6o and A Oq, the beaming direction A = corresponds to maximum polarization. The direction 
of maximum total intensity (i.e. the sum of intensities of both polarizations) does not coincide with the beaming direction, 
but even if we assume a detection threshold for a gamma ray polarimeter that is a function of the total intensity, the 
relative probability of detecting high polarization is considerably higher than in the "head-on" approximation. Writing 
the total intensity in the above limits, where relevant cp are at ^ <C 1, as 

/7r /'CO 
(1 -/3cos6l)-'=(l-|-cosw'^)smM'/>-- / (1 -/? cos 6')-^(l + coW^)^^, (17) 
-TT J —OO 

we find, performing the integrals, that 



For k=3, 



and, in a Euclidean space. 



J(A = 0) _,5/2 48 
/(A = i/r)-^ 83^^-^' 

VmaxiA = 0) ^ ^20) 



K.„,(A= 1/F) 

To summarize, scattering material in the Compton sailing state beams the polarized scattered component forward, 
and, for a roughly homogenous source distribution, a nearly fully polarized beam has the highest Vmax- The less intense 
unpolarized component, which by invariance of emitted power should be about as strong when averaged over a solid angle, 
presumably covers a larger solid angle, so for a source distribution with a low ^ , we might see mostly weakly-polarized 
GRB's. The above calculation has several simplifying assumptions: it is for the analytically tractable case 0oF » 1, 
A^ <C Go i.e. scattering off a hollow cone into a thin annulus. It neglects penetration to finite optical depths, multiple 
scattering, finite source size, and viscosity within the sheath. Many of these issues will be addressed in a subsequent paper 
with detailed numerical calculations (Levinson & Eichler astro-ph/0402457). 

Finally, it is worth recalling that the bulk Lorentz factor F of GRB fireballs, while usually estimated from observational 
constraints to be several hundred, has never been calculated from first principles. The phenomenon of Compton sailing 
suggests that a) there may be a connection between F and the geometry of the flow and b) F may vary continuously 
between the fireball and the ejected, supernova-associatc;d material from the envelope of the GRB's host star. It suggests 
the possibility of calculating F theoretically. While such a calculation would be formidable, the problem seems sufficiently 
well posed that F is no longer a mere free parameter imposed on the initial conditions of the fireball. 

I thank A. Levinson, E. Dcrishcv and Y. Lyubarsky for useful discussions. This research was supported by the Arnow 
Chair of Astrophysics at Ben Gurion University, by an Adler Fellowship awarded by the Israel Science Foundation, by 
a grant from the Israel-U.S. Binational Science Foundation and by a Center of Excellence grant from the Israeli Science 
Foundation. 
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